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Abstract 

The Wigner little group for massless particles is isomorphic to the Euclidean group SE(2). 
Applied to momentum eigenstates, or to infinite plane waves, the Euclidean “Wigner translations” 
act as the identity. We show that when applied to finite wavepackets the translation generators 
move the packet trajectory parallel to itself through a distance proportional to the particle’s helicity. 
We relate this effect to the Hall effect of light and to the Lorentz-frame dependence of the position 
of a massless spinning particle. 

PACS numbers: 03.30.+p, 03.50.De, 11.30.Cp 
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I. INTRODUCTION 


The Poincare group provides the fundamental kinematic symmetry of a relativistic parti¬ 
cle. As a non-compact group, all its unitary representations are infinite dimensional, but in 
a famous paper [T] Wigner showed that the physically interesting representations can be in¬ 
duced from finite-dimensional unitary representations of a little group , which is the subgroup 
of homogeneous Lorentz transformations that leaves some reference four-momentum in¬ 
variant. The representation space of the little group is the Hilbert space for the particle’s 
spin. If the particle has positive mass m, we may take as reference the four-momentum in the 
particle’s rest frame where p = (m, 0 ). The little group then consists of the space rotations 
SO(3). For a massless particle there is no rest frame and the reference momentum must be a 
null vector p—> (|p 0 |, po). The little group now consists of space rotations SO(2) about the 
three-vector po, together with operations that are generated by infinitesimal Lorentz boosts 
in directions perpendicular to po combined with compensating infinitesimal rotations. Re¬ 
markably the combined operations mutually commute, possess all the algebraic properties of 
Euclidean translations, and the resulting little group is isomorphic to the symmetry group 
SE(2) of the two-dimensional Euclidean plane. What is being moved by these translation 
operations? The answer given by Wigner is that they move nothing: if the translation gen¬ 
erators had a physical effect, the little-group representation would be infinite dimensional 
and the particle being described would have “continuous spin” — a property possessed by no 
known particle. Indeed the Wigner translations have no effect when applied to plane-wave 
solutions of the massless Dirac equation, and act as gauge transformations when applied 
to the vector potentials of plane-wave solutions of Maxwell’s equations [2j. Consequently 
they act as the identity on the momentum eigenstates created by the operator-valued coeffi¬ 
cients of the plane-wave modes, thus ensuring that the spin of a massless particle is entirely 
specified by a finite-dimensional representation of the SO(2) helicity subgroup [3j. 

It is the purpose of this paper to show that, while they have no effect on infinite plane 
waves, when applied to finite-size wave packets of non-zero helicity — and in particular to 
circularly polarized Gaussian packets — the Wigner translations do have an effect: they shift 
the wave packet trajectory parallel to itself. This shift is related to the relativistic Hall effect 
of light @H6] and to the observer dependence of the location of massless particles [7]. It gives 
rise to the unusual Lorentz covariance properties found mu in the chiral kinetic theory 
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approach to anomalous conservation laws pT0hTT2j and is also the source of the difficulty of 
obtaining a conventionally covariant classical mechanics for a massless spinning particle in 
a gravitational held mm- 

In section [XT] we will provide a suggestive algebraic argument for a sideways shift. In 


section III we will show that the shift actually occurs in finite-width beam solutions to 


Maxwell’s equations. In section IV we will discuss and resolve a potential paradox implied 
by the trajectory displacement. 


II. POINCARE ALGEBRA AND MASSLESS PARTICLES 


As an indication that Wigner translations can have a physical effect, we briefly review 
a well-known m H6] realization of the Poincare algebra for massless particles of helicity 
A in terms of quantum mechanical position and momentum operators. We start from the 
familiar commutators 

[xi,Pj\ = ihSij, \pi,Pj\= 0, (1) 

and use the fact that the spin of a massless particle is slaved to its direction of motion to 
motivate the definition of the angular momentum operator as 

Jk ^klm^lPm “1“ A . 'i • (2) 

IpI 

This unconventional definition preserves the usual commutation relation 

[A, A] ^fitklmPm- (3) 


However, in order to recover 


and 


\Jki A] 


( 4 ) 


[Jki Jl\ ( k! m Jm ■ 

we need to modify the commutator of the position-operator components to 


[Xk,Xl] 


■ f. a Pm 

iTiX &klm I ^ I q 

IpI 3 


( 5 ) 

( 6 ) 


Accepting that the position-operator components no longer commute, we can still use p° = 
|p| to define a generator of Lorentz boosts in direction k as 

K k = |(x fe |p| + |p|x fc ). (7) 
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These generators satisfy the remaining relations of the Lorentz Lie algebra 

^l\ ( 8 ) 

and act as expected on the momentum components: 

[K k) |p|] = ihp k , 

[K k ,Pi\ = ihS k i\p\. (9) 

We have therefore constructed a representation of the Poincare algebra on a quantum- 
mechanical ffilbert space. 

When we extend the algebra to include the position operators, things become more 
complicated. We find (at t — 0) 

[K k ,xi] = + |^| + ^ efcZm j^p| • ( 10 ) 

Neither term is immediately familiar. The expression in parentheses arises because the 
underlying Hamiltonian formalism automatically maintains the non-Lorentz invariant con¬ 
dition X 0 — t na. The term containing the hclicity A will be more interesting. 

We select a reference four-momentum p$ = (|Po|, Po) where p 0 = (0,0, p) and obtain 
the corresponding Wigner translation generators as the boosts and compensating rotations 
given by 


Ifi — K\ + J 2 , 

n 2 = A 2 — J\. (11) 


From pi) we see that these generators obey the SE(2) Lie algebra 


[ni,n 2 ] - o, [J 3 ,n x ] — zhn 2 , [^ 3 ,n 2 ] — — ihiii. 


( 12 ) 


From (J9J) and (10) we also see that x 1} x 2 , and the SE(2) generators collectively leave 
invariant the eigenspace with eigenvalues p = (0,0, p) and any fixed X 3 . Acting within the 
particular invariant subspace with x 3 = 0, we find that 


[U k ,xi] =-ihe k i 3 ~, (k,l = 1,2). 

p 


( 13 ) 
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In (13) the Wigner “translations” seemingly effect a genuine infinitesimal translation of the 


X\, X 2 coordinates in the X 3 = 0 plane, and hence a translation of the particle trajectory 
x(£) = (x\,X2,t) parallel to itself. Is this apparent displacement merely an artifact of an 
unconventional representation of the Poincare algebra, or does it have something to do with 
physics? 

In the next section we will use solutions of Maxwell’s equations to illustrate that this 
sideways shift is not just a mathematical curiosity, but corresponds to what occurs in nature 
- the trajectory of a circularly polarized photon is observer-dependent and is translated 
parallel to itself by an infinitesimal Lorentz boost and aberration-compensating rotation. 

III. PARAXIAL MAXWELL BEAMS 


We wish to consider the action of boosts and rotations on a finite-size photon wavepacket. 
It will serve to consider their effect on finite-width laser beam in the paraxial approximation. 
We will use units in which fj 0 = eo = c = 1. 

The scalar paraxial wave equation 


d 2 x d 2 x , w .dx 

w 2 + WT = 0 

ox z oy z oz 


is obtained from the full scalar wave equation 

d 2 (j) d 2 (j) d 2 (j) d 2 (p 


+ 


+ 


dx 2 dy 2 dz 2 dt 2 


= 0 


(14) 


(15) 


by writing 


<j>(T,t)=x(r)e ik ^ (16) 

and assuming that y(r) is sufficiently slowly varying that we can ignore its second derivative 


d 2 x/dz 2 in comparison to the remaining terms in (J14[). 


The simplest solution of eq. (14) is the Gaussian-beam US] 


x(r) = 


] r I „,2 I „,2 


(z - iz 0 ) 


exp 


x +y~ , x + y 

+ ik- 


2 w 2 (z) 


2 R(z) 


where 


z 2 z 2 z 2 z 2 

wHz) = R(z) = 


(17) 


(18) 


kz 0 ' ' z 

In this solution the beam is propagating in the +z direction, the quantity w(z) is the width 
of the beam at a distance z away from its waist, and R(z) is the radius of curvature of 
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FIG. 1: Slice through a paraxial scalar beam with parameters k = 10, zq = 10. a) Density plot of 
original beam amplitude Re{x(®j 0, z, b)e ikz }; b) Beam amplitude after Lorentz transformation (eq. 


(26)) with rapidity s = 0.5; c) Beam amplitude after both Lorentz transformation and aberration- 
compensating rotation though 9 = — tan^ 1 (sinh s) = —31.5°. 


the wavefront passing through the point r = (0, 0, z). The width grows linearly with 2 once 
0 zq , and the angular half-width is l/kw(0). The condition for the paraxial approximation 
to be accurate (kzo 1) is equivalent to the beam having small asymptotic divergence. We 
will always be interested in the region z < Zq where the beam is narrow and almost parallel 
sided. 


From any two independent solutions /, g of the scalar paraxial equation we can find [III] 
vector E and B fields that are internally consistent solutions of Maxwell’s equations up to 
accuracy of order l/(kl) 2 , where l is some charateristic length such as zq 


and 


E x {r,t) 

E y (r,t) 

E z {r,t) 


1 fd 2 f 8 2 f\ 1 d 2 g 

i (r ’ + 4 k 2 V dx 2 dy 2 ) + 2 k 2 dxdy ’ 

1 fd 2 g d 2 g\ 1 d 2 f 

9 ' T,t ' ~ 4 k 2 1 ~dx? ~dy 2 ) + 2 k? dxdy’ 

i_(df dg\ 

k \dx dy) ’ 


B x {r,t) 

B y (r,t) 

B z (r,t) 


1 fd 2 g d 2 g\ 1 d 2 f 
+ 4 Ik? V^2 ~d^) + 2k? dxdy ’ 

ft D 1 ( d2 f 92 1 92 9 

/lr ’ ’ 4 k 2 \dx 2 dy 2 ) 2k 2 dxdy’ 

_i(dg_ df\ 

k \dx dy J 


(19) 


( 20 ) 
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To obtain a Gaussian TEMqo beam that is circularly polarized with positive helicity we take 


/(r, £) = y(r)e* fc ( 2_t ) and g(r,t) = ix(r)e lk( ^ z ~ t \ with y(r) given by eq. (17). 

Using Mathematica™ to manipulate the resulting rather lengthy expressions we find, for 
example, that the time-average energy density in the beam is 


T 00 = -(|E | 2 + | B | 2 ) = ^ + ^ + ^ + Z D) 2 c -kz a (x 2 W)l(z 2 +zl) 


8(^ 2 + £q ) 3 


( 21 ) 


and the three components of the time-averaged Poynting vector S = (E x B) are 

c _ T io _ 0 3 ^ - x2 v z o + x v 2z - y 3z o + 4 ( xz - v z q)(t 2 + 4)) -fc* 0 (* a +y a )/(* a +*g ) 

1 2(^ 2 + ^ 0 2 )3 

q _ T 20 _ ( y 3 z + y 2xz 0 + X 2 yz + x 3 z 0 + 4 (yz + xz 0 )(z 2 + zg)) bzo(x 2 W)/(z 2 +z 2 n) 

y 2(z 2 + z 2 ) 3 

q _ T so _ (■~ x4 - 2x2 y 2 - y A + 16 (- 2 + ^o) 2 ) ^- fen ( 3; 2 + ^ )/(z 2 +2 2 ) 
z 8(^ 2 + ^ 2 ) 3 
The energy-flux streamlines twist in the direction of the beam helicity [2U], consequently the 

z component of the angular momentum density 


( 22 ) 


Yj z xSy yS x 


(23) 


is non-zero. If we integrate over the plane z = 0we find that 


and 


p <U f 

1 Z - 


j def 
J z = 


z=0 


n(-l + 8k 2 ZQ ) 2 

S ’ dxdy = —W4 — = kZ , <1 + 0 


E ,dxdy = 


2 = 0 


7 r(l + 2kzo) 

k 3 z% 


-\zr« {l+0 


( kzo ) 2 

1 


(24) 


(25) 


.(^o) 2 . 

The ratio P z /J z is equal to k in region (. kzo 1) where paraxial aproximation is accurate. 
This is what is to be expected: P z gives the linear momentum per unit length, which should 
be hk per photon; J z gives the angular momentum per unit length of the beam, which should 
be h per photon. 

We now compute the E and B fields as seen from a reference frame moving along the +x 
axis at rapidity s. The corresponding Lorentz transformation takes 


E x (x,y,z,t ) E x (x',y,z,t'), 

E y (x,y,z,t) H > E y (x', y, z, t') coshs — B z (x', y, z, t', ) sinhs 
E z (x,y, z,t) e-x E z (x', y, z, t') coshs + B y (x\ y, z, t') sinhs, 
B x (x,y,z,t ) B x (x',y,z,t') 

B y (x,y, z,t ) H > B y (x', y, z, t') coshs + E z (x', y, z, t') sinhs, 
B z (x,y,z,t ) i —> B z (x'y, z, t') coshs — E y (x', y, z,t') sinhs, 
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(26) 











where 


x — x cosh s + t sinh s 

t' — t cosh s — x sinh s. (27) 

The Lorentz transformation changes the wave vector from k = (0,0, k) to k' = 
(k sinh s', 0, k), so the direction of propagation has been rotated though an aberration angle 
of \6\ = tan -1 (sinh s). The wavefronts are therefore tilted. The beam envelope, however, 
still lies parallel to the £-axis, and is moving towards the observer at speed j3 = tanh s (see 
figure [l}b). 

The Lorentz transformation also affects the energy density distribution and the Poynting- 
vector flux though the z = 0 plane. In addition to a Lorentz contraction it noticeably shifts 
the position of their maxima (see figure [2j) . To quantify these shifts we can compute the 
location of the Lorentz transformed energy density and energy flux centroids. The required 
integrals are still Gaussian and can be done analytically. With the definition 


£ = 


T 00 dxdy, 


(28) 


2=0 


we have 


[Ay] density — ^ 


2=0 


yT w dxdy , 
zo(4 + 8 kzo sinlis) 


(1 + 8 kzo + 8k 2 Zq) cosh s — 4 kzo sech s ’ 


and 


= L anh T 1 + o \(fc 0 y. 


[Ay]fl ux = [[ yS z dxdy, 

*z JJz= 0 

2 ^ 0 (1 — 2kzo) tanh s 
1 — 8k 2 zl 


(29) 


= — tanh s <( 1 + O 
2 k 


(30) 


(fe o) 2 . 

For positive hclicity, both centroids are displaced to the left when seen from an observer 
moving towards the upward-propagating beam. The centroids do not coincide, the energy- 
flux centroid moving only half as far as the energy-density centroid. Such displacements 
are not restricted to Gaussian beams. A similar boost-induced sideways shift and centroid 
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separation was exhibited in [7] for Bessel beams possessing orbital angular momentum. It 
was also explained there that the centroid separation arises solely from the geometrical 
effect pointed out in [2T]: because of their corkscrew trajectories, energy-flux streamlines 
passing through a surface rotated away from perpendicular to the direction of propagation 
find themselves inclined at different angles to the surface to the right and left of the plane 
of rotation. Consequently, even in the absence of a Lorentz boost, the energy-flux centroid 
of a tilted beam is displaced with respect to its energy-density centroid [ 21 j. 

We wish to obtain a finite-displacement version of the Wigner translations, so, after 
performing the boost, we rotate the Lorentz transformed beam about r = 0 though an 
aberration-compensating angle of tan _1 (sinh .s). After the rotation the wavevector becomes 
k = (0,0, k coshs) and the wavefronts again lie parallel to the x-y plane. Consequently 
the energy-flux streamlines no longer possess a left-right asymmetry. We find numerically 
that the position of the energy centroid in the z = 0 plane is unchanged by the rotation 
(T 00 is a scalar under space rotations) while the energy-flux centroid moves into coincidence 
with the energy-density centroid. Thus, as result of the combined boost and compensating 
rotation both centroids have been shifted though a distance Ay = (1 /k) tanh s = (3/ k, where 
f3 = v/c. The beam spot is restored to its pre-boost appearance, and we could repeat the 
operation and translate the beam spot through a further distance. If we reverse the hclicity, 
we change the sign of this shift. 

In the absence of the lateral shift, the combination of boost and compensating rotation 
would leave the trajectory of a short wavepacket emitted from r = 0 at t = 0 unchanged. The 
continuous beam, which can be though of as arising from a stream of sequentially emitted 
wavepackets, is not left invariant, however. How it changes is shown in fig. [l]-c. We see that 
the transformed beam can be though of as a sequence of pulses each fired in the +z direction 
by an emitter that is moving rapidly to the left. It is reminiscent of a diagonal steam of 
strictly upward-moving projectiles fired from a horizontally moving gun in the old Atari™ 
game “Space Invaders.” Any particular packet continues to move parallel to the 0 axis, but 
as a result of the lateral shift in the z = 0 plane, its entire trajectory is shifted sideways by 
Ay = (l/fc)tanhs. Figure [l]-c also shows why the action of the Wigner translations take 


their simple form (13) only in the plane £3 = z = 0. In any other plane the translations get 
mixed up with the geometric effect of the rotation. 

The finite-s boosts considered in this section have effects on the photon energy and 
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FIG. 2: Beam spot profiles in the z = 0 plane for k = 10, zq = 3. a) Original intensity T 00 (x,y); 


b) Lorentz transformation ofT 00 (x,y ) under eq. (26) with rapidity s = 2.0. The spot center is at 
y = 0.095; c) Poynting energy flux S z (x, y) = T 30 after Lorentz transformation. The spot center is 
at y = 0.0475; d) Poynting energy flux S z (x. y) = T 30 (x, y) after aberration-compensating rotation. 
The flux maximum is at y = 0.095. The rotated intensity distribution has similar appearance, and 
its maximum is also at y = 0.095. 


intensity that appear at quadratic order in the rapidity s. If we alternate a sequence of 
infinitesimal boosts and compensating rotations, the quadratic terms can be neglected and 
only the sideways shift (now equal to A /p times the total rapidity change) remains. We are 
in effect assembling a Trotter-product approximation that converges to exponentials of the 
Wigner translation generators 
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z c) 



FIG. 3: a) A pair of massless particles with spin S collide head-on; b) The particles viewed from 
above in a frame moving towards the collision; c) A front view from the moving frame shows the 
particles miss one another. 

IV. DISCUSSION 

The direction and magnitude of the boost-induced lateral shift can be understood from 
a geometric picture (See [8] for a related argument). Consider two massless particles, both 
possessing helicity p-S spin /|p| = A and heading directly towards one another parallel to the x 
axis. Because they will collide head-on, they have no relative orbital angular momentum and 
the two spin angular momenta S sp i n = (±A, 0, 0) also sum to zero. Seen from a frame moving 
along the y axis towards the collision point, however, the unit vectors in the direction of 
the particles’ motion have components (isechs, — tanhs, 0). Because the spin of a massless 
particle is slaved to its direction of motion there is now a net spin component of 2A tanhs 
directed towards the observer. Nonetheless, in the new frame, the total angular momentum 
will remain zero so the spin contribution must be offset by an orbital angular momentum 
of —2Atanhs. This orbital angular momentum can only come from a lateral shift of each 
particle’s trajectory by Ac = (±A/|p|) tanhs (see fig. [3]). For a photon p = h k and A = h, 
so we recover the shift seen in our Gaussian beam. Of course, if two particles collide and 
produce two pions in one frame they must produce two pions when seen from another 
frame. That the particles apparently miss each other because of the sideways shift cannot 
affect the pion production. The incipient paradox is resolved by the fact that partial-wave 
scattering amplitudes depend only on the total relative angular momentum J = L + S sp i n 
of the particles, and this quantity is not affected by the shift. The shift still has physical 
consequences, though. If we move a detector such as a photographic emulsion though the 
beam, it will be sensitive to either the energy density or the energy flux in its own rest 
frame, and these quantities have been displaced by the motion. 


12 












V. ACKNOWLEDGEMENTS 


This work was supported by the National Science Foundation under grant number NSF 
DMR 13-06011. In the course of this work MS has exchanged many useful emails with 
Peter Horvathy and Christian Duval. MS would also like to thank Konstantin Bliokh for 
discussions about the observer dependence of the location of light rays, and for drawing our 
attention to [Tj. We also thank the authors of j8] for sending us an early version of their 
work. 


[1] E. Wigner, On unitary representation of the inhomogeneous Lorentz group, Ann. Math. 40, 
149-204 (1939). 

[2] D. Han, Y. S. Kim, D. Son, Gauge transformations as Lorentz-boosted rotations, Physics 
Letters 131B, 327-329 (1983); Photons, neutrinos, and gauge transformations, Am. J. Phys. 
bf 54, 818-821 (1986). 

[3] S. Weinberg, Feynman Rules for Any Spin. II. Massless Particles, Phys. Rev. 134, B882-B896 
(1964); Photons and Gravitons in S-Matrix Theory: Derivation of Charge Conservation and 
Equality of Gravitational and Inertial Mass, 135, B1049-B1056 (1964); See also The Quantum 
Theory of Fields, Vol. 1, Cambridge University Press, Massachusetts (1996). 

[4] M. Onoda, S. Murakami, N. Nagaosa, Hall Effect of Light, Phys. Rev. Lett. 93, 083901 (2004). 

[5] C. Duval Z. Horvath, P. A. Horvathy, Geometrical SpinOptics and the Optical Hall Effect, 
Journ. Georn. Phys. 57, 925 (2007). math-ph/0509031. 

[ 6 ] C. Duval, Z. Horvath, P. A. Horvathy, Fermat Principle for polarized light, Phys. Rev. D74, 
021701 (R) (2006). cond-mat/0509636. 

[7] K. Y. Bliokh, F. Nori, Relativistic Hall effect, Phys. Rev. Lett. 108, 120403 1-4 (2012). 

[ 8 ] J-Y. Chen, D. T. Son, M. A. Stephanov, H-U. Yee, Y. Yin, Lorentz Invariance in Chiral 
Kinetic Theory, Phys. Rev. Lett. 113, 182302 (2014). 

[9] C. Duval, M. Elbistan, P. A. Horvathy, P.-M. Zhang Wigner-Souriau translations and Lorentz 
symmetry of chiral fermions, ai'Xiv:1411.6541 

[10] M. A. Stephanov, Y. Yin, Chiral Kinetic Theory, Phys. Rev. Lett. 109, 162001 1-5 (2012). 

[11] M. Stone, V. Dwivedi, A Classical Version of the Non-Abelian Gauge Anomaly, Phys. Rev. 


13 






D88, 045012 1-8 (2013). 

[12] V. Dwivedi, M. Stone, Classical chiral kinetic theory and anomalies in even space-time di¬ 
mensions, J. Phys. A: Math. Theor. 47, 025401 1-20 (2014). 

[13] M. Stone, V. Dwivedi, T. Zhou, Berry Phase, Lorentz Covariance, and Anomalous Velocity 
for Dirac and Weyl Particles , arXiv: 1406.0354 

[14] C. Duval, P. Horvathy, Chiral fermions as massless spinning particles, arXiv:1406.0718, 

[15] M. V. Atre, A. P. Balachandran, T. R. Govindarajan, Massless spinning particles in all di¬ 
mensions and novel magnetic monopoles, Int. J. Modern Physics A 2, 453-483 (1987). 

[16] B.-S. Skagerstam, Localization of Massless Spinning Particles and the Berry Phase, hep- 
th/9210054, 

[17] A. J. Hansson, T. Regge, C. Teitelboim, Constrained Hamiltonian Systems, (Accademia 
Nazionale dei Lincei, Roma 1976) p33. 

[18] A. Garg, Classical Electromagnetism in a Nutshell, Princeton University Press (2012) pl64. 

[19] W. L. Erikson, Surendra Singh, Polarization properties of Gaussian laser beams, Phys. Rev. 
E49, 5778-5786 (1994). 

[20] M.V.Berry, K.T McDonald, Exact and geometrical optics energy trajectories in twisted beams, 
J. Opt. A: Pure Appl. Opt. 10, 035005 1-7 (2008). 

[21] A. Aiello, N. Lindlein, C. Marquardt, G. Leuchs, Transverse Angular Momentum and Geo¬ 
metric Spin Hall Effect of Light, Phys. Rev. Lett. 103, 100401 1-4 (2009). 


14 





